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The magnetic ordering of superconducting single crystals of K0.85Fe1.83Se2.09 has been studied 
between 10 K and 550 K using ^'^Fe Mossbauer spectroscopy. Despite being superconducting below 
Tsc ~ 30 K, the iron sublattice in Ko.85Fei,83Se2.o9 clearly exhibits magnetic order from well below 
Tsc to its Neel temperature of Tjv = 532 ± 2 K. The iron moments are ordered perpendicular to the 
single crystal plates, i.e. parallel to the crystal c-axis. The order collapses rapidly above 500 K and 
the accompanying growth of a paramagnetic component suggests that the magnetic transition may 
be first order, which may explain the unusual temperature dependence reported in recent neutron 
diffraction studies. 

PACS numbers: 



I. INTRODUCTION 

The co-existence of magnetism and superconductivity 
in the new iron-chalcogenide superconductors raises the 
possibility of unconventional pairing mechanisms that 
may be associated with their magnetism. [H-Q As with 
the more established cuprate superconductors, the iron- 
based superconductors have layered structures; the pla- 
nar Fe layers tetrahedrally coordinated by As or chalco- 
gen anions (Se or Te) are believed to be responsible for 
superconductivity. Stacking of the FeAs building blocks 
with alkali, alkaline earth or rare earth-oxygen spacer 
layers forms the basic classes of iron arsenic supercon- 
ductors in these compounds: Ill-type AFeAs [5|, 122- 
type AFeaAsafi-li, 1111-type ROFeAs [13, [H and more 
complex block containing phases, e.g. Sr2V03FeAs [l^ . 
SrgScaFesAsaOg [11, Sr4Sc2Fe2As206 [3- The simple 
binary 11-type iron chalcogenides have no spacer layers 
and superconductivity can be induced by doping FeTe 
with S[15j] or Se[l5. Unlike the other iron-based super- 
conductors, FeSe is a superconductor [l3|, Tsc 8 K, 
with no static magnetic order and its transition temper- 
ature can be increased up to 37 K by applying pressure 
[Tsf or 15 K in FeSeo.sTeo.s [H. More recently, supercon- 
ductivity above 30 K has been reported in AxFe2_ySe2 
(A = K, Cs, Rb or Tl) [MH, a compound with the 
same unit cell structure as the AFe2As2 compounds. 

Recent neutron diffraction measurements on 
K0.8Fe1.6Se2 have suggested that not only does magnetic 
order co-exist with superconductivity but that the iron 
moments are remarkably large (3.31 fis/Fe) and their 
magnetic ordering temperature may be as high as 560 K 
[2^ . The relatively complex antiferromagnetic structure 
places all of the iron moments parallel to the c-axis. 
Muon spin relaxation (pSR) has confirmed the high 
magnetic ordering temperature in Cso.8(FeSeo.98)2 (25[ 
but the development of a paramagnetic component near 
Tn suggests that the transition may be first order in 



nature rather than being a more conventional second 
order transition. First order magnetic transitions are 
commonly associated with changes in crystal structure, 
and both synchrotron x-ray diffraction [26j and neutron 
diffraction [13, H^l have now shown evidence for a 
structural change from 14/ m to M/mmm associated 
with a disordering of iron vacancies that occurs in 
the vicinity of the magnetic transition. However, the 
magnetic structure adopted by the iron sublattice has 
not been confirmed, and room temperature neutron 
diffraction studies of CsyFe2-xSe2 [1^ and AyFe2_xSe2 
(A — Rb, K) '28^ have suggested that the iron moments 
may be much smaller (~2.5 /is/Fe) and also that the 
magnetic structure may be far more complex than 
initially suggested, with the iron atoms being distributed 
among two (magnetically) inequivalent sublattices and 
carrying very different magnetic moments. Moreover, 
even the ordering direction has been questioned and it is 
possible that the iron moments may lie in the ab-plane, 
at least for CsyFe2_xSe2 |26|, ra ther than parallel to the 
c-axis as initially suggested [2J] . 

Given the many questions surrounding the magnetic 
ordering of the iron moments in the AyFe2_xSe2 sys- 
tem, we have undertaken a ^^Fe Mossbauer study 
K0.85Fe1.83Se2.09- While Mossbauer spectroscopy cannot 
be used to determine magnetic structures directly, it is 
a quantitative local probe that can be used to set hard 
limits on possible structures. As we will show below, 
the observation of a single, well-split magnetic compo- 
nent allows us to rule out any structure in which the 
iron sub-lattice is further subdivided into multiple, in- 
equivalent sites, and the scale of the splitting is incon- 
sistent with smaller moment estimates. Furthermore, 
by working with single crystal samples, we are able to 
demonstrate that the moments order parallel to the c- 
axis, ruling out any models that invoke planar ordering. 
Finally, the development of a paramagnetic component 
that co-exists with the magnetically ordered phase and 
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that grows at its expense on heating through T^r, con- 
firms that the magnetic transition is indeed first order in 
nature. 



II. EXPERIMENTAL METHODS 

The preparation and characterisation of the crystals 
used in this study has been described in detail elsewhere 
[23I I so only a basic description will be provided here. 
Single crystals of K0.85Fe1.83Se2.09 were grown from a 
Ko.8Fe2Se2 melt, as described in Ref . [20| . First the FeSe 
precursor was prepared by reacting stoichiometric Fe and 
Se at 1050°C. Then, K and FeSe with a nominal compo- 
sition of Ko.8Fe2Se2 were placed in an alumina crucible 
that was sealed in an amorphous silica tube. The growth 
was placed in a furnace in a vented enclosure and heated 
to 1050°C, where it was held for a 2 hour soak. The fur- 
nace temperature was then slowly lowered to 750°C over 
50 hours; the furnace was then turned off and the sample 
"furnace cooled" over an additional 10 hours. Once the 
ampoules were opened, large (~ 1 x 1 x 0.02 cm^) dark 
shiny crystals could be mechanically separated from the 
solidified melt. 

Crystals were characterised by powder x-ray diffrac- 
tion using a Rigaku Minifiex X-ray diffractometer. The 
average chemical composition was determined by exam- 
ining multiple points on a cleaved surface of the crystal, 
using wavelength dispersive x-ray spectroscopy (WDS) in 
a JEOL JXA-8200 electron microscope. A backscattered 
electron analysis (BSE) was performed using an acceler- 
ating voltage of 20 kV. Magnetic susceptibility was mea- 
sured in a Quantum Design MPMS, SQUID magnetome- 
ter. The x-ray diffraction pattern can be indexed using 
the I4/mmm space group, with lattice parameters refined 
by Rietica of a = 3.8897(8)A and c = 14.141(3)1, in 
good agreement with previous values [13]. The in-plane 
resistivity of the furnace cooled sample is very similar to 
that of earlier reports (2^ 2M ■ ^ broad resistive maxi- 
mum centered near 160 K is followed by a drop of nearly 
a factor of 6 {psooK / P35k)- The sharp transition to zero 
resistance gives an onset temperature of Tsc = 30.1 K. 
Corresponding features were also seen in the susceptibil- 
ity and heat capacity. 

The as-grown crystals were too thick to serve as 
Mossbauer absorber but we found that they could be 
cleaved quite easily using a razor blade, much like mica, 
and many large-area plates could be formed from each 
crystal. Two single crystal mosaic samples were pre- 
pared from the same batch of crystals. The first, for 
low-temperature work, was prepared by attaching sev- 
eral single crystal plates to a 12 mm diameter disc of 
100/im thick Kapton foil using Apiezon N grease. Care 
was taken to ensure that there were no gaps, but minimal 
overlap between the crystals. This sample was trans- 
ferred promptly to a vibration-isolated closed-cycle re- 
frigerator with the sample held in vacuum. The second 
sample, for the high-temperature work, was attached to a 



2-inch diameter 10-mil beryllium disc using diluted GE- 
7031 varnish before being mounted in a resistively heated 
oven, again with the sample in vacuum. While we oper- 
ated somewhat above the maximum service temperature 
of the varnish, the sample was cycled above 250°C three 
times without any evidence of degradation. 

The Mossbauer spectra were collected on conventional 
spectrometers using 50 mCi ^^CoRh sources mounted 
on electromechanical drives operated in constant accel- 
eration mode (on the high-temperature system) and sine- 
mode (on the low-temperature system). The spectrom- 
eters were calibrated against a— Fe metal at room tem- 
perature. The closed-cycle refrigerator cools to 10 K, 
with temperature sensing and control using a calibrated 
silicon diode mounted on the copper sample stage. Mea- 
sured gradients (centre to edge of sample) in the oven 
are less than 1 K up to 750 K. Control and sensing rely 
on four type-K thermocouples. Temperature stability in 
both cases is better than 0.2 K. Spectra were fitted using 
a conventional non-linear least-squares minimisation rou- 
tine to a sum of equal-width Lorentzian lines. Magnetic 
patterns were fitted assuming first-order perturbation in 
order to combine the effects of the magnetic hyperfine 
field (Bhf) and the electric field gradient. 

III. RESULTS 

Several conclusions can be reached simply by inspec- 
tion of the spectrum taken at 10 K (Fig. [Ij. The spec- 
trum is dominated by a single, well-split, magnetic com- 
ponent. This confirms that K0.85Fe1.83Se2.09 is indeed 
magnetically ordered in the superconducting state (re- 
call Tsc ^ 30 K for this sample). The single magnetic 
component allows us to rule out any magnetic structures 
involving multiple iron sub-sites with moments that dif- 
fer by more than a few percent. As we will show be- 
low, the large hyperfine field (Bhf ^29 T) is inconsistent 
with a small iron moment and so places further limits 
on possible magnetic structures. Finally, two of the lines 
normally present in a magnetically split ^^Fe Mossbauer 
spectra, are clearly absent from the 10 K pattern. 

A magnetic field at the ^^Fe nucleus, either externally 
applied or transferred from an ordered moment on the 
iron atom, lifts the degeneracy of the nuclear states and, 
in combination with the selection rules for the f — > 5 
transition, leads to a six-line pattern with intensities of 
3:R:1:1:R:3 (counting from left to right in Fig. [1]). For a 
powder sample, R=2, however if there is a unique angle, 
0, between the magnetic field and the direction of the 
7— beam used to record the spectrum, then the intensity, 
R, of the Ami = transitions is given by: 

4 sin^ e 

1 + cos2 e 

R=0 implies that 6 is also zero so that the magnetic field, 
and by extension, the moments that lead to it, is paral- 
lel to the 7— beam. Since the sample consists of an ab- 
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FIG. 1: (color online) ^^Fe Mossbauer spectra of 
K0.85Fe1.g3Se2.09 showing the evolution of the magnetic or- 
dering on heating from 10 K (well below Tsc ~ 30 K) to 
533 K where the material is paramagnetic. The absence of 
the Ami transitions in the ordered state indicates that the 
moments are parallel to the crystal c-axis (see text), while the 
growth of a central paramagnetic component above 500 K is 
characteristic of a first order magnetic transition. Solid lines 
are fits as described in the text. 



plane mosaic of single crystals, setting 6 — means that 
the magnetic ordering direction is parallel to the c-axis, 
ruling out any magnetic structures that involve planar 
ordering of the iron moments. We note that R is a rel- 
atively soft function of 9 near zero, and a free fit to the 
intensity of the Am/ — transitions is consistent with an 
angle of up to 20°, but this does not lead to a significant 
improvement in for the fit, nor would such an angle be 
consistent with a planar ordering of the iron moments. 
Indeed, if the ordering were planar, then R would be 4, 
and the Am/ = transitions would provide the strongest 
features in the spectrum. 

Estimating the iron moment from the observed hyper- 
fine field requires some care as the scaling is imperfect 
at best [sot. However, some data exist on binary iron- 
chalcogenides that can be used as a guide (Table HI . If 
we use the factor of 6.2 T/fiB for FctScs with our mea- 



sured Bhf of 29.4 T we obtain a rather large estimate of 
4.7 /is/Fe for the iron moment in this system. This is 
significantly larger than the 3.31 /xs/Fe reported on the 
basis of neutron diffraction '24], however it does suggest 
that the iron moment is indeed substantial as even the 
larger conversion factor for the sulphide yields 3.5 /is/Fe. 
If we assume that B^/ is at least proportional to the iron 
moment, then we can use the observed change in Jihf be- 
tween 10 K and 293 K to scale the 3.31 ^s/Fe observed 
at 11 K[23| to get an estimate of 3.0 /is/Fe for the mo- 
ment at room temperature for comparison with the much 
smaller value of 2.55 /is/Fe reported by Pomjakushin et 
al.^^. However, the strong temperature dependence of 
magnetic signal noted by Bao et al.^^ suggests a very 
rapid decline in ordered moment to about 2.8 iis/Fe by 
room temperature. It is possible that much of the vari- 
ation may be intrinsic to the material and its variable 
stoichiometry, so that combined measurements on a well 
characterised sample will be needed to settle this. 

TABLE I: Average hyperfine fields (B^/) derived from *^Fe 
Mossbauer spectroscopy and average iron moments derived 
from neutron diffraction for approximately equi-atomic iron- 
chalcogenide compounds with estimated field-moment con- 
version factors. The Fe-Te system exhibits significant vari- 
ability and measurements have yet to be made on com- 
mon samples making the conversion factor unreliable. There 
is however a clear trend to lower values in the sequence 
S^Se-^Te. 



Compound 


Average 
(T) 


Average 
moment 
Ms/Fe 


Conversion 
Factor 


Sulphides 
FerSs 


26.8f32] 


3.16[33] 


8.5 


Selenides 
FerSeg 


24.1 [34] 


3.86 [35] 


6.2 


Tellurides 
Fei.iasTe 
Fei+xTe 

0.076 < X < 0.141 

Fei.oesTe 

Fei.osTe 

Fei.iiTe 

Fei.osTe 


11[40| 
10.34f41] 


2.07[36] 
1.96-2.03[37] 

2.25[38] 
2.54[39] 


4.3-5.2 



Impurities may provide a possible origin for the vari- 
ation in measured moments. Mossbauer spectroscopy, 
while sensitive to the presence of impurity phases, does 
not rely on normalisation to the total sample in order to 
determine moments, they come rather from the observed 
line splitting, and not the intensity. Neutron diffraction, 
by contrast, while providing far more information on the 
magnetic ordering, ultimately relies on peak intensities, 
normalised to the total nuclear scattering, to determine 
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the magnetic moments. It is clear from the 10 K spec- 
trum shown in Fig. [I] that there is a central paramagnetic 
component present that involves about 12±2% of the iron 
in the sample. Such high apparent impurity levels in sin- 
gle crystal samples with no impurities detected by pow- 
der x-ray diffraction[23[, deserves further attention. If 
the paramagnetic component is not an "impurity" then 
it must either be intrinsic to the structure or a property 
of the material. 

At the temperatures of interest here, K0.85Fe1.g3Se2.09 
adopts a vacancy-ordered lA/m modification of the par- 
ent ThCr2Si2— type lA/mmm structure with iron essen- 
tially filling a 16z site and lea ving ordered vacancies on 
the (almost) empty 4d site [1^ l3ll | Occupations of '^8% 
for the Ye-Ad site have been reported [Sri. If we assume 
full occupation of the Fe-16i site in our sample, this leaves 
12.5% of the iron in the Ad site, while scaling the com- 
position to two selenium atoms per formula unit gives a 
lower estimate of 9% for the iron in the Ad site. Partial 
occupation of the Fe-16i site would leave more iron to be 
accommodated in the Ad site. As we see no evidence for 
a second magnetic component that could be associated 
with iron in the Ad site, it is possible that the iron in 
these more isolated sites does not order, in which case 
our estimate of 9-12.5% in the Ad site is fully consistent 
with the 12±2% paramagnetic component observed in 
the Mossbauer spectrum. 




whether this spatial variation is only associated with the 
surface or is representative of the bulk behaviour of the 
sample. It should be noted, though, that such patterns 
appear in samples grown by furnace cooling as well as 
samples decanted from a liquid melt [2^. Regardless of 
its origin, the non-magnetic Fe signal represents a mi- 
nority of the Fe sites and does not substantially change 
with temperature until the first order phase transition is 
reached. 
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FIG. 3: (color online) Temperature dependence of the mag- 
netic hyperfine field (Bh/) in K0.85Fe1.83Se2.09. The solid line 
is a fit to a J=| Brillouin function between 200 K and 500 K 
that yields an expected transition of 600±30 K, well above 
the observed value of 532±2 K. Fitted errors on Yihf a-re less 
than 0.1 T, much smaller than the plotting symbols. The 
rapid collapse above 500 K is accompanied by the growth of 
a paramagnetic component (see Fig. [4]). 




FIG. 2: Backscattered electron analysis (BSE) image of a 
cleaved crystal surface of K0.85Fe1.83Se2.09 taken at an ac- 
celerating voltage of 20 kV. The lighter regions have lower 
potassium concentrations than the darker background area. 

Another possible origin of the 12±2% non-magnetic Fe 
component in the low temperature (including room tem- 
perature) state can be seen in the backscattered electron 
analysis (BSE) image shown in Fig. [2] This image re- 
veals that there is, at the micron scale, a modulation in 
the surface composition that can be correlated, through 
a preliminary line-scan analysis of the WDS data, with 
reductions of K content in the lighter regions. Given the 
length scale associated with these regions, combined with 
the probing volume of the WDS analysis, it is not clear 



Raising the temperature leads to the expected decline 
in B?i/, however it is clear from Fig.[T]that magnetic order 
persists up to 530 K, confirming that K0.85Fe1.83Se2.09 
has a remarkably high ordering temperature. The tem- 
perature dependence of B/j/ shown in Fig. [3] yields an 
ordering temperature of Tat = 532 ± 2 K. However this 
is not the result of the fit to a J=^ Brillouin function 
shown in Fig. |3] as this predicts a transition temperature 
of 600±30 K and the observed behaviour departs from 
this curve above 500 K. The two points that bracket the 
transition are at 530 K, where a clear magnetic signal is 
seen, and at 533 K where the sample is no longer mag- 
netic, setting the transition at 532 ± 2 K. 

A neutron diffraction study of K0.8Fe1.6Se2 found two 
regions in which the temperature dependence of the mag- 
netic parameter was unusual [IJ. From 50 K to 450 K 
they found a linear dependence of the (101) magnetic 
peak intensity, suggesting that is a linear function 
of temperature. The clear curvature of B/j/(T) in this 
region, shown in Fig. [31 is not consistent with this form, 
as squaring our observed Bft/(T) to get something that 
would scale with the scattering intensity in a neutron 
diffraction pattern leads to increased curvature rather 
than linear behaviour. 

Above 500 K, Wei Bao et al. reported a very rapid 
decrease in the (101) intensity [IJ] leading to an order- 
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FIG. 4: (color online) Temperature dependence of the mag- 
netic fraction in K0.85Fe1.83Se2.09- The rapid collapse above 
500 K indicates that the magnetic transition has first order 
character and may be associated with a structural transition. 



ing temperature of '-^^560 K. While our sample compo- 
sition is slightly different and our ordering temperature 
slightly lower, we see the same abrupt loss of magnetic 
order in Fig. [H Inspection of the spectra above 500 K 
shown in Fig. [1] reveals that the intensity of the mag- 
netic peaks decreases visibly as their splitting falls. The 
ability to uniquely separate the amount of a magnetic 
phase (seen through line intensities) from the magnitude 
of the magnetic order (seen independently through line 
splittings) is an important strength of Mossbauer spec- 
troscopy. Tracking the fraction of the iron that is present 
as a magnetically ordered form (Fig. 2]) confirms that 
the magnetic phase is disappearing even faster than the 
splitting that marks the order. This strongly suggests 
that the magnetic phase is transforming before it reaches 
its true ordering temperature (which we estimate to be 
about 600 K) and that the observed transition is being 
driven by a first order structural event. This view is 
supported by the neutron diffraction work of Wei Bao et 
al. [2JI where they also tracked the intensity of the (110) 
structural peak that is associated with the /4/to vacancy- 
ordered structure of K0.8Fe1.6Se2 below 580 K. This peak 
starts to lose intensity at the same temperature at which 



the (101) magnetic peak starts its sudden decline. As we 
see both a weakening of the magnetic order and a reduc- 
tion in the magnetic fraction above 500 K, it is possible 
that the break-up of the vacancy-ordered magnetic form 
reduces the magnetic connectivity of the ordered phase 
until it forms a non-percolating network of finite clusters. 
The magnetic order is then lost at a temperature below 
both its intrinsic ordering temperature, and the temper- 
ature at which the vacancy-ordered 14/ m structure fully 
transforms to the high-temperature lA/mmm form. 
IV. CONCLUSIONS 

Our ^""Fe Mossbauer spectroscopy study of single crys- 
tals of Ko.ssFei. 83862. 09 confirms the presence of mag- 
netic order from well below Tsc ^ 30 K to Tjv = 
532 ± 2 K. The large magnetic splitting of 29.4±0.1 T at 
10 K indicates that the iron moments are large, consistent 
with values of 3.31 /Fe observed by neutron diffraction 
at 11 k[23|, while the line intensities indicate that the 
ordering is parallel to the c-axis. An apparent paramag- 
netic impurity phase is attributed to iron atoms in the 4c? 
site. Analysis of the spectra taken in the vicinity of Tjv 
shows that the magnetic fraction decreases rapidly above 
500 K and that the loss of order is driven by a first order 
structural transition as the material transforms into the 
high-temperature lA/mmm form. 



Acknowledgements 

Financial support for various stages of this work was 
provided by the Natural Sciences and Engineering Re- 
search Council of Canada and Fonds Quebecois de la 
Recherche sur la Nature et les Technologies. JMC ac- 
knowledges support from the Canada Research Chairs 
programme. R.H. and P.C.C. are supported by AFOSR- 
MURI grant #FA9550-09- 1-0603. W.E.S., S.L.B. and 
P.C.C. are also supported by the U.S. Department of En- 
ergy, Office of Basic Energy Science, Division of Materials 
Sciences and Engineering. Ames Laboratory is operated 
for the U.S. Department of Energy by Iowa State Uni- 
versity under Contract No. DE-AC02-07CH11358. 



[1] Kenji Ishida, Yusuke Nakai, and Hideo Hosono, J. Phys. 
Soc. Jpn., 78 062001 (2009). 

[2] M. D. Lumsden and A. D. Christianson, J. Phys.: Con- 
dens. Matter, 22 203203 (2010). 

[3] Paul C. Canfield and Sergey L. Bud'ko, Annual Review 
of Condensed Matter Physics, 1, 27 (2010). 

[4] Johnpierre Paglione and Richard L. Greene, Nature 
Physics 6, 645 (2010). 

[5] X. C. Wang, Q. Q. Liu, Y. X. Lv, W. B. Gao, L. X. Yang, 
R. C. Yu, F. Y. Li, C. Q. Jin, Solid State Commun. 148, 
538 (2008). 



[6] M. Rotter, M. Tegel, and D. Johrendt, Phys. Rev. Lett. 
101, 107006 (2008). 

[7] Athena S. Sefat, Rongying Jin, Michael A. McCuire, 
Brian C. Sales, David J. Singh, and David Mandrus, 
Phys. Rev. Lett. 101, 117004 (2008). 

[8] N. Ni, M. E. Tillman, J.-Q. Yan, A. Kracher, S. T. Han- 
nahs, S. L. Bud'ko, and P. C. Canfield, Phys. Rev. B 78, 
214515 (2008). 

[9] A. Leithe-Jasper, W. Schnelle, C. Geibel, and H. Rosner, 
Phys. Rev. Lett. 101, 207004 (2008). 
[10] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, 



6 



J. Am. Chem. Soc. 130, 3296 (2008). 
[11] X. H. Chen, T. Wu, G. Wu, R. H. Liu, H. Chen, and D. 

F. Fang, Nature 453, 761 (2008). 
[12] X. Zhu, F. Han, G. Mu, P. Cheng, B. Shen, B. Zeng, and 

H. H. Wen, Phys. Rev. B 79, 220512 (2009). 
[13] X. Zhu, F. Han, G. Mu, B. Zeng, P. Cheng, B. Shen, H. 

H. Wen, Phys. Rev. B 79 024516 (2009). 
[14] G. F. Chen, T. L. Xia, H. X. Yang, J. Q. Li, P. Zheng, 

J. L. Luo, N. L. Wang, Supercond. Sci. Tech. 22, 072001 

(2009). 

[15] Rongwei Hu, Emil S. Bozin, J. B. Warren, and C. Petro- 

vic, Phys. Rev. B 80, 214514 (2009). 
[16] Y. Mizuguchi, Y. Hara, K. Deguchi, S. Tsuda, T. Yam- 

aguchi, K. Takeda, H. Kotegawa, H. Tou and, Y. Takano, 

Supercond. Sci. Tech., 23 054013 (2010). 
[17] F. C. Hsu, J. Y. Luo, K. W. The, T. K. Chen, T. W. 

Huang, P. M. Wu, Y. C. Lee, Y. L. Huang, Y. Y. Chu, 

D. C. Yan and M. K. Wu, Proc. Nat. Acad. Sci. 105, 

14262 (2008). 

[18] S. Medvedev, T. M. McQueen, L Trojan, T. Palasyuk, M. 
L Eremets, R. J. Cava, S. Naghavi, F. Casper, V. Kseno- 
fontov, G. Wortmann and C. Felser, Nature Mater., 8 
630 (2009). 

[19] Jiangang Guo, Shifeng Jin, Gang Wang, Shunchong 

Wang, Kaixing Zhu, Tingling Zhou, Meng He and Xi- 

aolong Chen, Phys. Rev. B 82, 180520 (2010). 
[20] J. J. Ying, X. F. Wang, X. G. Luo, A. F. Wang, M. 

Zhang, Y. J. Yan, Z. J. X iang, R. H. Liu, P. Cheng, G. 

J. Ye, X. H. Chen. larXi^ :1012.5552vl (2010). 
[21] Chun-Hong Li, Bing Shen, Fei Han, Xiyu Zhu, Hai-Hu 

Wen, arXiv:1012.5637v2 (2010). 
[22] Minghu Fang, Hangdong Wang, Chiheng Dong, Zu- 

juan Li, Chunmu Feng, Jian Chen, H. Q. Yuan, 

larXiv:1012 .5236 (2010). 
[23] R. Hu, K. Cho, H. Kim, H. Hodovanets, W. E. 

Straszheim, M. A. Tanatar, R. Prozorov, S. L. Bud'ko, 

P. C. Canfield, arXiv:1102.1931 (2011). 
[24] Wei Bao, Q. Huang, G. F. Chen, M. A. Green, D. M. 

Wang, J. B. He, X. Q. Wang, Y. Qiu, larXiv:1102.0830l 

(2011). 

[25] Z. Shermadini, A. Krzton-Maziopa, M. Bendele, R. 

Khasanov, H. Luetkens, K. Conder, E. Pomjakushina, 

S. Weyeneth, V. Pomjakushin, O. Bossen, A. Amato, 

larXiv:ll0l .l873vl (2011). 
[26] V. Yu. Pomjakushin, D. V. Sheptyakov, E. V. 

Pomjakushina, A. Krzton-Maziopa, K. Conder, D. 

Chernyshov, V. Svitlyk, Z. Shermadini larXiv: 1101. 19191 

(2011). 

[27] Wei Bao, G. N. Li, Q. Huang, G. F. Chen, J. B. He, M. A. 
Green, Y. Qiu, D. M. Wang, J. L. Luo, arXiv:1101.3674 
(2011). 

[28] V. Yu. Pomjakushin, E. V. Pomjakushina, A. Krzton- 
Maziopa, K. Conder, Z. Shermadini, I arXiv: 1101.33801 
(2011). 

[29] D. M. Wang, J. B. He, T.-L. Xia, G. F. Chen, 

arXiv:1101.0789vl (2011). 
[30] S.M. Dubiel, J. Alloys Compounds 488, 18 (2010). 
[31] P. Zavalij, W. Bao, X.F. Wang, J.J. Ying, X.H. Chen, 

D.M. Wang, J.B. He, X.Q. Wang, G.F Chen, P-Y Hsieh, 

Q. Huang, M.A. Green, arXiv:1101.4882 (2011). 
[32] H. Kobayashi, M. Sato, T. Kamimura, M. Sakai, H. On- 

odera, N. Kuroda and Y. Yamaguchi, J. Phys.: Condens. 

Matter, 9, 515 (1997). 
[33] A.V. Powell, P. Vaqueiro, K.S. Knight, L.C. Chapon, and 



R.D. Sanchez, Phys. Rev. B 70, 014415 (2004). 
[34] H.N. Ok and S.W. Lee, Phys. Rev. B 8, 4267 (1973). 
[35] A.F. Andresen and J. Leciejewicz, J. de Phys. 25, 574 

(1964). 

[36] D. Fruchart, P. Convert, P. Wolfers, R. Madar, J.P. Sen- 

ateur and R. Fruchart, Mater. Res. Bull. 10, 169 (1975). 
[37] Wei Bao, Y. Qiu, Q. Huang, M.A. Green, P. Zajdel, M.R. 

Fitzsimmons, M. Zhernenkov, S. Chang, M. Fang, B. 

Qian, E.K. Vehstedt, J. Yang, H.M. Pham, L. Spinu, 

and Z.Q. Mao, Phys. Rev. Lett. 102, 247001 (2009). 
[38] S. Li, C. de la Cruz, Q. Huang, Y. Chen, J.W. Lynn, J. 

Hu, Y.-L. Huang, F.-C. Hsu, K.-W. Yeh, M.-K. Wu, and 

P. Dai, Phys. Rev. B 79, 054503 (2009). 
[39] A. Martinelli, A. Palenzona, M. Tropeano, C. Ferdeghini, 

M. Putti, M.R. Cimberle, T.D. Nguyen, M. Affronte, and 

C. Ritter, Phys. Rev. B 81, 094115 (2010). 
[40] E. Hermon, W.B. Muir, J. Quaroni, and R.C. Sweet, 

Can. J. Phys. 52, 1800 (1974). 
[41] Y. Mizuguchi, T. Furubayashi, K. Deguchi, S. Tsuda, T. 

Yamaguchi, Y. Takano, Physica C 470, S338 (2010). 



